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The high precision attained by cosmological data in the last jew years has 
increased the interest in exact solutions. Analytic expressions for solutions 
in the Standard Model are presented here for all combinations of A = , 
A / 0, k = and k ^ 0, in the presence and absence of radiation and 
nonrelativistic matter. The most complete case (here called the K^CDM 
Model) has A ^ 0,/t ^ 0, and supposes the presence of radiation and dust. 
It exhibits clearly the recent onset of acceleration. The treatment includes 
particular models of interest such as the ACDM Model (which includes the 
cosmological constant plus cold dark matter as source constituents). 
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1 INTRODUCTION 



Recent cosmological data [U El E] suggest a present-day Universe dominated 
by dark energy, though with a significant contribution of matter (visible and 
dark). This brings to the fore a Friedmann- Robertson- Walker model with 
a strong de Sitter component, that is, with a cosmological constant. In the 
usual approach to standard cosmology the history of the Universe is divided 
into eras, each one dominated by one of the source-constituents. Separate 
solutions are used for each era. The intermediate periods were described by 
approximations which were quite acceptable until recently. The precision of 
present-day data, however, makes it desirable to have solutions as complete 
and exact as possible. An exact solution for the Standard Model can in prin- 
ciple be obtained provided an equation of state (giving the pressure in terms 
of the energy density) is given for each constituent. For all times except 
the first few seconds (the equation of state for a very-very-high-temperature 
black-body is unknown), the main difficulty concerns dark matter. It is our 
aim here to provide a short, compact presentation, in terms of present-day 
parameters, of the solutions [H El [6] for a dust dark matter, which can be 
called "KryCDM Model" . This complete, unified version is a necessity when- 
ever applications and comparisons are in view. We shall, however, proceed 
step by step, the reason for that being that the general, complete solution 
is actually implicit (gives time as a function of the scale parameter) and 
involves elliptic functions of the first and the third kinds. Such a progres- 
sive presentation seems advisable, describing solutions for particular values 
of the cosmological constant and curvature parameters, presence and absence 
of matter and radiation, with increasing complexity. 

Section El is a short introduction to the Friedmann equations [3 El El E3ZI1 
[TTj . the main objective being the introduction of convenient parameters and 
notation. Maybe too much detail on well-known results is given in section El 
concerned with de Sitter solutions. The assumptions there made concerning 
inflation will, however, propagate to the ensuing solutions. In effect, we use 
the standard exponential inflation there found as a kind of "correspondence 
principle": all other solutions will generalize that case, and reduce to it for 
the corresponding choice of parameters. Section H] treats the model with 
radiation, cosmological constant and curvature. The limits to particular 
cases, as the pure radiation model, are readily obtained. In section El we 
analyze the case for matter plus another contribution, either the cosmological 
constant or curvature. Adding any other contribution leads to a solution as 
involved as the general case, which is given in Section El 
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2 GENERAL OUTLOOK 



Let us start fixing notation by recalling Einstein's equations with a cosmological- 
constant and source energy content modeled by a perfect fluid: 

8vrG 



R 



T 

„4 ^' 



T, 



(p + pc 2 ) u^Uy - p 



Here p = e/c 2 is the mass equivalent of the energy density. We use metric sig- 
nature (+, — , — , — ), so that UjjU^ 1 = 1 for timelike flux lines. Spacetimes with 
homogeneous and isotropic space sections are described by the Robertson- 
Walker line element [3, E] 

dr 2 



ds 2 = c 2 dt 2 - a 2 (t) 



+ r 2 d6 2 + r l sin 2 



(2) 



1 — KT 1 

with k = 0, ±1, which reduces Einstein's equations to the two Friedmann 
equations for the scale parameter a(t): 
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Combining the last two equations, one finds the expression for energy 
conservation, 

dp 3 ( p 
da 

The total matter energy density e 
relativistic matter ("baryons"), dark or not, and radiation. It is convenient 
to relate to A the "dark energy" density 

Ac 4 

eA = ^ ■ (6) 



p m c 2 = PbC 2 + p-yC 2 includes both non- 
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In terms of the Hubble function 
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Eqs. (l3"P| assume the forms 
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Actually, these equations acquire their most convenient form in terms of di- 
mensionless variables, obtained by dividing by the present-day value H 2 {t = 
to) = Hq and introducing the parameters 



n K {t) = - 
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with pcrit = P^?. The first equation, in particular, becomes 
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which leads to the usual normalization for present-day values, 
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Observation will fix the value of k: we are forced to choose k = +1 if Qtotai > 1 
and k = — 1 if Qtotai < 1- I n both cases the present-day value of the scale 
factor is determined by ao = <J Ut t K ; _ 1 "^- If Qtotai = 1, then k = and ao 

remains undetermined. 

Pressure in equation (jHJ) is the sum p = p m = Pb + P-y, the necessary 
expressions being given by thermodynamic considerations. Non-relativistic 
matter is represented as an ideal gas, i.e., pb = nkT. This pressure, however, 
is negligible against its related density (dust approximation): (p b +Pb/c 2 ) = 
^(mc 2 + kT) ~ pb, since mc 2 ^> kT. The radiation equation of state is that 

of a black body 
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The conservation equation (J5J), which is 



dpb g Pfc 
da a 



+ 



dp 1 + 4 P 1 
da a 



0, 



(15) 



can be solved immediately under the hypothesis of independence of the com- 
ponents: matter and radiation have their own separated dynamics with re- 
spect to the scale factor a. The result is the pair 



Pb 



Pob = Pob (1 + zf ; p 7 = p 07 (^j = po 7 (1 + zf (16) 
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where z is the red-shift and po — p{ a o)- These results for the relationship 
between density and scale factor for matter and radiation is a natural con- 
sequence of the baryonic number conservation, which is strongly supported 
by the very high value of the proton lifetime [13]. In fact, if the number of 
baryons N b = n b a 3 is set constant, -r^ = 0, it follows that n b a 3 = n^al, n b 
being the numerical baryonic density. The first of equations (1161) then follows 
because p b = n b (m b c 2 ), and leads to the second once inserted into (IT51) . 
Finally, we introduce the notations 

2 8ttG 4 2 4 f-\r,\ 

7C = -^-po y a = n j0 H a , (17) 

Mc 2 = ^p 0b a 3 = n b0 H 2 al , (18) 
r 2 Ac 2 

^ = ^ = ^o 2 , (19) 
in terms of which equations ( |3|4l) take on the simple forms we shall be using: 

a 2 c 2 7c 2 Mc 2 kc 2 

a = c 2 __7c^_Mc 2 _ 

a L 2 a 4 2 a 3 ' 1 J 

In proceeding to our step-by-step approach, it will be convenient to start 
from the simplest and dominant solution, de Sitter spacetimes. They will help 
fixing some aspects to be retained in the other cases. Notice that ( l20l I2TT) 
impose severe conditions on parameters for the particular solutions we shall 
be examining. An example which we shall bypass is the case A = 7 = M = 0, 
which would only be admissible for k — 0, — 1. 



3 DE SITTER SOLUTIONS 

We recall that there are two different kinds of such spaces: 

1. de Sitter spacetime proper: indicated dS: a one-sheeted hyperbolic 4- 
space with topology R 1 x S 3 , an inclusion in the pseudo-Euclidean 
space E 4 ' 1 satisfying (A, B = 0, 1, 2, 3, 4) 

Vab Z A Z B = Vap Ce ~ (e 4 ) 2 = - L 2 (22) 
Its scalar curvature R is — within our conventions — negative. 
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2. anti-de Sitter spacetime: indicated AdS: a two-sheeted hyperbolic 4- 
space with topology S 1 x R 3 , an inclusion in the pseudo-Euclidean 
space E 3,2 satisfying 

It has R > 0. With the notation 7/44 = s, dS and AdS can be put together in 
flABt A Z B = fhfi + = *L* . (24) 

Both dS (s = -1,A > 0) and AdS (s = +1,A < 0) are solutions of the 
sourceless Einstein's equations with a cosmological constant related to the 
pseudo-radius L and the scalar curvature by 

A = "4 = -f. (25) 

As they are sourceless solutions and have homogeneous and isotropic space 
sections, these spacetimes satisfy equations ( !20l) and ( l2~TT) with vanishing 7 
and M: 

/ r \ 2 



a 2 = - s \ -J a 2 - kc 2 ; (26) 

a = -5 (^) 2 a • (27) 

Let us start with = 0. The second equation has solutions of the 

general form a(t) = A e at + B Q e^, actually with a = f3. Thus, the general 
expression reduces to the form a(t) = A e at . But then a/a = H = a = 
±^/— s f ■ The only way to have a real positive constant H is to choose s < 
and the upper sign. These choices lead to 

a{t)=Ae s t\ (28) 

the standard inflationary solution of exponential type. The initial condition 
is set for the beginning of the time measure, i.e., A = a(t = 0). We could 
equivalently choose as "initial" value the present time to, such that A = 
aoe~ 



Notice that A = — s-p->0is essential for inflation — only the dS space 
(and not the AdS) can lead to inflation. Another point: the usual treatment 
considers sometimes solutions with initial value A = a(0) = 0. Here, however, 
there will be no inflation if the initial value of a(t) is zero. We shall from 
now on consider only A > 0, A > 0, and look for solutions generalizing this 
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result, that is, which reduce to flzHj) when only a cosmological term is present 
and k — 0. 

For k^O, the solution allowing for inflation is 



a(t) 



A cosh ~ + vM 2 - kL 2 sinh f 
1 
2 



A + VA 2 - kL 2 ) e l + \A - VA 2 - kL 2 I 

If we introduce the function 

f{t) = [A + VA 2 -kL 2 ^) e"i 



and use the relation 



A - VA 2 - kL- 
solution ff29l) takes on the form 

a(t) = 

Summing up, the solution is 



kL 2 



(A + VA 2 - kL 2 ) 



(k^O) 



(29) 



(30) 



(31) 



f(t) 



A = ^^0: a(t) = 

L 2 [f(t)=(A+ VA 2 - kL 2 ) ei 



kL 2 
W) 



m + % 



(32) 



with 

ct 



(33) 



and the Friedmann-Robertson- Walker form of the inflationary de Sitter line 
element will consequently be 



ds z = c A dt 



/(*) 



kL 2 
W) 



dr 2 
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+ r 2 (d6 2 + sin 2 9d4> 2 



(34) 



4 RADIATION WITH k ^ 0, A ^ 



From all the evidence we have today, the Universe has passed an initial stage 
during which, besides the cosmological term, ultrarelativistic matter or radi- 
ation provided a significant contribution as a source in Einstein's equations. 
This era is described by equations f[2"01 and fl2"TT) with M = 0, 
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(35) 
(36) 
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We solve the first of these equations by direct integration (completing 
squares and changing variables) with our conventional initial condition a(0) = 
A. The result reads: 



a(t) 



kL 2 
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1 
2 



g(t) + 



git) 



(37) 



with the definition 



9(t) 



kL 2 



A 2 ) - V7L 2 + A 2 (A 2 - kL 2 ) 



2 — t 



(31 



The solution above involves, besides the cosmological constant and k, 
other quantities hidden in parameter 7: the present-day value ao of the 
expansion parameter and the actual value of the radiation density p 7 o- At 
t = we recover the initial condition a(0) = A, as it should be. 

It is a good test of consistency to analyze some limits, (i) In the absence 
of radiation, 7 = 0, solution f[3"Tj) reduces dutifully to the de Sitter solution for 
k^O, equation ( !29l) . (ii) In the limit of no cosmological constant, L — > +00, 
the scale factor ( 1371) is simplified to 



a(t) = y A 2 + 2^7 - kA 2 ct - K(ct) 2 . 
(iii) Taking only k = in (1371) . we are left with: 



(39) 



a(t) 



-| A 2 + v / 7^TA4 e 2 r*+ ^-7^+^ e - 2 r*|. (40) 



(iv) If we, in addiction the restriction k = 0, perform the limit A — ► (i.e. 
L — > 00), it follows 

(41) 



a{t) = ^/A 2 + 2^ct, 

the conventional solution for a radiation dominated universe, in which a oc 
t 1 / 2 . 



5 NON-RELATIVISTIC MATTER 

This section deals with the possible cosmological models in absence of radi- 
ation: 7 = 0. It is divided in two special cases: first, besides M/0, one 
assumes A ^0 but k = 0; and, then, one considers no cosmological constant, 
A = 0, but k ^ 0. 



S 



5.1 Matter with A^O,k = (ACDM Model) 

Equation (jHJ) can be used to eliminate the term e m from 
approximation p = leads then, for k — 0, to 



3 /c 2 
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2 V^ 2 

whose integration is straightforward: 

[i + H(0)] e^ T [i-H(0)] 



The dust 



(42) 



L[f + H(0)] e^±li-H(V 



(43) 



The integration constant H (0) is obtained from (1201) in terms of previously 
defined measurable quantities: 



# 2 (0) 



Mc 2 J_ 
~A^ + L 2 



The expansion parameter follows from expression (120]) for H 



(44) 



a(t) = (Mc 2 ) 



or, substituting ( |43|) with the choice of the upper sign, 
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e h' 



(45) 



The choice of sign is adequate, as it leads to the appropriate limits: (i) 
M -> reduces (HS) to (j2g|); and (ii) with L -> oo (i.e. A -> 0) we find 



2 

a (t) = (a 3 / 2 + ?V^c*Y , 

the solution for a matter dominated universe. In terms of still another con- 
venient function, 

3 c 



J(t) 



ML 2 
A 3 



- — t 



(46) 



solution ( |45l) becomes 



a(t) = A j(i) 



ML 2 1 
Wj^t) 



(47) 
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5.2 Matter with A = 0, k ^ 

Equations ( !20l) and ( I2T1) are, in this case, 



, n MC 2 

a 2 + kc 2 = 

a 

Mc 2 



(48) 
(49) 



The first equation, 



ada 



\J Ma — na 1 



c l dt 



(50) 



is better solved separately for each value of k: 

(I) Case k = can be used as a test: integration gives a 3 ' 2 — 

A 3/2 



Met, just the result of the previous section. 

(II) Case k > 0, actually k — 1: a change of variable puts the integrals 
into the form 



a ^/ydy 



1 1 y_ 

x M 



Mc I dt. 

o 



(51) 



Integration gives a(t) in an implicit way: 



ct = VMA - A 2 - VMa - a 2 + 

M 



arcsm 




M 



— arcsm 




(52) 



(III) Case k < 0, actually k = — 1: integration now gives 



ct = y/Ma + a 2 - VMA + A 2 + Mm 



(53) 



When M > 0, a function of ^ turns up: 



>\/m + ( m ) 



£L _|_ /I _i_ _£L 



M + V 1 + M 



(54) 
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This can be rewritten as 



7-1/ \ P V^Vl+ 

F(x) 



y/x + yJl + X 

(55) 



aft) 
Af 



4 1 ct 

1 e« 



Function F(x) has actually a very simple, monotonic behavior, and so has 
a(t). 

We have in this section solved Friedmann's equations with a dust source 
(M 0), choosing solutions of inflationary type. Actually, only the cases 
with either A ^ or k ^ have been presented. The case with both A ^ 
and k is much more involved. In fact, all the remaining classes of 
solutions are practically as involved as the general case including all the four 
components A, k, M and 7. Let us then proceed directly to that case. 

6 GENERAL CASE: 7 ^0,M^0,A^0,^0 

This A-fCDM case is described by the complete equations (|20f2ip . We repeat 
one of then: 

a[t) +KC a(t) 2 a(t) L 2 " U - (5bj 
We shall integrate this first-order equation, which is equivalent to 



a(t) d(t) = \j C 1^L + 7C 2 + Mc 2 a(t) - Kc 2 a(t) 2 (57) 



or 

c r* . f a a da 



dt= / . (58) 
L Jo J A ^fa A + 7L 2 + ML 2 a - kL 2 o? 

The procedure to arrive at the solution is long and rather cumbersome. 
Let us proceed step by step: 

1. Rewrite the integral above in terms of the roots {rj} of the denominator: 

f a da f a da 



y/a* + 7L 2 + ML 2 a - KL 2 a 2 J y/(a - n)(a - r 2 )(a - r 3 )(a - r 4 ) 
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2. Introduce some constants, in terms of which the roots will be expressed 
later on: 

W = [27(ML 2 ) 2 -2(kL 2 ) 3 + 72 7 kL 4 ] 2 , 



U = W+ ^W 2 -A[{kL 2 ) 2 + l2 1 L 2 f 
V = l{ 2 «L*-(|) 5 -(iy [(.L 2 ) 2 + 12 7 L 2 ] 
X = (2kL 2 -V). (59) 

3. The roots will then be: 

r 3 = -\ (W+^FT^) , r 4 = -| . (60) 

4. Solutions will be implicit, and involve elliptic functions of two kinds [T4"l 

da us]: 

• the elliptic integral of first kind with parameter m and amplitude 

0: 

m] = /" 1 d0. (61) 

•/ o yl - m sin 

• the elliptic integral of third kind with parameter m, characteristic 
n and amplitude 0: 

II[0, n, m] = [ — - M. (62) 

Jo (1 — nsin 9)yl — msin 6* 

5. The characteristic and the parameter turning up will be: 

n. (63) 



/ r 2 


- r 4 \ 




- r 3 \ 


n = \ — 




; m = 




v n 


- r 4 y 


\ r 2 


- r 3 J 



6. The amplitude turning up will be 



fa) = arcsin \ l —. ■ (64) 

n(a — t\) 
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7. The implicit solution will then be 



ct (7*1 — r 2 ) {Tl[0(a), n, m] — II [0(A), n, m]} 

2L \/( r 2 ~ r 3)( r i - r 4) 

r 1 {F[<f ) (a),m]-F[<p(A),m]} 

\Z(r 2 - r 3 )(ri - r 4 ) 

We see from their definitions that both elliptic functions vanish when 
= 0. The intricacy of this solution justifies its use only in the general case. 
It does reduce to the expected particular cases: this can be more easily ver- 
ified by confronting the plots of each particular case with the corresponding 
graphics constructed from the general scale factor obtained from (R?5|) . 



Model 
A y CDM 

Radiation 

Matter 

A CDM 



t 

2.5 5 7.5 10 12.5 15 17.5 

Figure 1: The scale factor a(t)/ao in terms oft (in Gyear). The two higher 
curves are reference models (pure radiation and matter, respectively) exhibit- 
ing an overall deceleration. The full curve comes from the A'yCDM universe, 
which coincides with the ACDM model in the scale used. It presents an in- 
flection point, at which A starts dominating and acceleration sets on. 




The graphs shown in Figure [T] include the general solution (full line, the 
A'yCDM Model), and two other reference cases: the radiation dominated 
Universe (dashed line) and the matter dominated Universe (dotted curve). 
Following the observational data, the curvature parameter is set to zero for 
all functions a(t), k = 0. The values of L, M and 7 are chosen in accordance 
with the present day values P [13]: ft to tai = 1-02 ± 0.02, Q w = 0.27 ± 0.04 
(contribution of baryonic and dark matter) and f2 7 o = (4.9 ± 0.5) x 10~ 5 , 
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and Hq = (71I00Q3) km sec -1 Mpc -1 . Time T — ct (c = 1) is given in 
Gyear and the scale factor is measured in units of its present day value a , 
an arbitrary constant. When considering pure radiation (pure matter), we 

take f2 7 o = fitotal (^60 — ^total)- 

The A'yCDM model gives 13.8 Gyear for the age of the universe. Unlike 
the two other models, it exhibits an inflection point for not too large a z, 
which starts the present day acceleration. As expected, the effect of matter 
and radiation surpasses largely that of the cosmological constant in the early 
Universe, and a pure radiation model evolves much faster than a matter 
dominated universe. In the present and future, however, the cosmological 
constant becomes more important than radiation and matter and dominates 
the cosmic dynamics. 

Another interesting plot is that of matter and A (with 7 = k — 0), the 
ACDM model. It is in fact also plotted in Figure ED, but coincides with 
the complete solution in the scale used. At small scales, however, it is clear 
that the general A^CDM model {Vt m = 0.27, fi A = 0.73, 7 = k = 0) differs 
from the ACDM model, which corresponds to the lower curve in Figure [21 
which is an enlargement of Figure [T] for early times. The radiation is the 
dominant component in this primeval period: see the dashed curve close to 
the vertical axis, or the interval (t < 10 kyear) where complete model (full 
line) overcomes the matter one (dotted line). 

7 FINAL COMMENTS 

The solutions given above hold for all times, as long as dark matter is a dust 
formed by baryons or, at least, when its energy density behaves according to 
the pure volumetric law e m = e m0 (1 + z) 3 . Notice that Eq. (l2"Tj) is, in time 
units of Hq 1 , 

- = n A - n j0 (1 + zf - \ n w (i + zf, (66) 

(X 

so that acceleration will change sign when 

Ov = 0,0 (1 + zf + \ n b0 (1 + zf . (67) 

This condition is k- independent. For the present-day favored [3] values Qa = 
0.73, f2&o — 0.27, f2 7 o = 4.88 x 10~ 5 , this take-over happens rather recently, 
at z ~ 0.75, or a(t) ~ 0.57 a®. The limitations of the A'-fCDM model can 
be seen in the estimate value coming up for the cross-over, which seems too 
high: (1 + z) = Oo/Oyo ~ 5530, to be compared with the preferred WMAP 
value z w 3234. 
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Figure 2: Zooming of the curve a(t)/ao for the primeval universe, with time t 
now given in kyear. The dashed curve of the radiation model lies very close 
the vertical axis. The difference between '-/ACDM and ACDM curves (full 
line and lower dashed line) is clearly shown. 



We have seen that the ACDM model, described by equation ( l45l) . has a 
simpler form than the general A'yCDM solution, equation ( 1651) : the latter 
gives the scale factor as an implicit function of time, unlike the former. Figure 
[U shows also that the ACDM model coincides with the complete solution at 
recent times. Nevertheless, the curve of A'yCDM model differs from the one 
of ACDM in the early universe (Figure [2]), when radiation has a dominant 
role. Whenever precision is needed, the general solution (I65D is to be used. 
This general solution will be valid as long as we take the combination of ultra- 
relativistic and non-relativistic ideal particles in the presence of a vacuum 
energy with equation of state p = — e. 
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